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Summary 


Methods for relating or mapping estimates of volumetric bone mineral density (VBMD) obtained via 
quantitative computed tomography (QCT) to bone volume fraction (BVF) are outlined mathematically. 
The methods are based on definitions of bone properties, cited experimental studies, and regression 
relations derived from those studies for trabecular bone in the proximal femur. BVF values calculated 
from four different methods were compared with the experimental average and numerical range of values 
obtained from the intertrochanteric region of male and female human subjects, aged 18 to 49. The BVF 
values computed from the conversion methods used data from two sources. One source provided pre-bed- 
rest VBMD values in the intertrochanteric region from 24 bed rest subjects who participated in a 70-day 
study. Another source provided preflight VBMD values from 18 astronauts who spent 4 to 6 months on 
the International Space Station (ISS). To aid the use of a mapping from BMD to BVF, the discussion 
includes how to formulate the conversions for the purpose of computational modeling. An application of 
the conversions would be used to aid in computational modeling of time-varying changes in VBMD as 
they relate to changes in BVF via bone remodeling and/or modeling. 


Introduction 


Most computational models of bone remodeling or bone adaptation track changes in apparent density 
(Refs. 1 to 3) or bone volume fraction (BVF) (Refs. 4 to 6), but in vivo measurements of bone density are 
given in terms of bone mineral density (BMD). Areal BMD (aBMD) is defined as bone mineral content 
divided by total area, g/cm”. Volumetric BMD (vBMD) is bone mineral content divided by total volume, 
g/cm>. Dual-energy x-ray absorptiometry, an enhanced form of x-ray technology, is an established 
standard bone density scanning method for measuring aBMD (Ref. 7). One medical technique that 
measures VBMD is quantitative computed tomography (QCT) (Ref. 8). QCT uses a standard x-ray 
computed tomography (CT) scanner with a calibration standard to convert Hounsfield units of the image 
to VBMD values. In order to compare model prediction to experimental data, a method of converting 
aBMD or vBMD to BVF is needed. 

In NASA’s Digital Astronaut Project (DAP), the musculoskeletal research component has as a current 
task the development of a computational bone physiology model for space flight bone physiology 
analysis (Ref. 9). One aim is to understand the changes in vVBMD, particularly bone loss, at various bone 
sites and the response to in-flight and postflight exercise countermeasures. The mathematical formulation 
describes the removal and replacement of bone via bone remodeling and simulates changes in bone in 
terms of changes in BVF. The changes are observed in subjects of bed rest studies that simulate 
microgravity and in astronauts during and after 4- to 6-month periods on the International Space Station 
(ISS). In addition, the DAP bone physiology model can add higher fidelity by breaking up BVF into 
mineralized volume fraction, M, and osteoid volume fraction, O: 
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BVF =M+O (1) 


Because of its three-dimensional capability and its ability to provide specific data for the separate cortical 
and trabecular regions of bone, the DAP bone physiology model is developed directly for QCT data. 
Human QCT scans are generally lower resolution to minimize radiation health risks, which in turn tends 
to prevent the determination of BVF via image processing software. Therefore, a method for converting 
initial VBMD data to initial BVF is needed in order to simulate time-course changes in BVF that can then 
be converted back to VBMD. This study presents a number of conversion methods with a focus on the 
trabecular region of the proximal femur. Terms used in this paper are defined in the Appendix. 


Acronyms 
aBMD areal bone mineral density 
BMD bone mineral density 
BVF bone volume fraction 
CFT countermeasures and functional testing 
CT computed tomography 
DAP Digital Astronaut Project 
ISS International Space Station 
mCT microcomputed tomography 
QCT quantitative computed tomography 
vBMD volumetric bone mineral density 
Symbols 

aBMD areal bone mineral density 
BMD bone mineral density 
BVF bone volume fraction 
M mineralized volume fraction 
O osteoid volume fraction 
vBMD volumetric bone mineral density 
or ash fraction o = Pash 

Papp 
p density 


Relations and Formulas 


Algorithms for the conversions are based in part on using linear correlation relations obtained from 
experimental studies that relate QCT-derived VBMD, pqcr , to ash density, Pash, and ash density to 
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apparent dry density, Papp. From experimental results on proximal femur samples, conversion from pact 


tO Pash for combined cortical and trabecular developed by Keyak et al. (Refs. 10 to 11) is given by 


Pash = 0.887pQct + 0.0633 (2) 


The experiments did not include results of correlating Papp to Pash prior to ashing for the proximal 
femur. However, an earlier work by Keyak et al. (Ref. 12) involving the proximal tibia trabecular region 
produced the relation 


Papp = 1-66 Pash + 0.00457 (3) 


Similar experimental work by Schileo et al. (Ref. 13) using human femur specimens obtained the 
following regressions for trabecular specimens (4a), which was not significantly different from the 
regression line for cortical specimens (4b): 


Papp =1-64Pash + 0.01 (4a) 


Papp =1-58Pash + 9.11 (4b) 


In addition to the experimentally derived regressions shown above, key experimental studies have 
been used to obtain relational estimates of true density and BVF. Experimental work from Keller 
(Ref. 14) computed Pash, Papp, and ash fraction a from 199 spine specimens and 297 femur specimens. Ash 
fraction values ranged from 0.174 to 0.662 with a mean spine value of 0.610 and a mean femur value of 
0.583. Hernandez et al. (Ref. 15) concluded from that study that the ash fraction of unmineralized osteoid, 
a = 0, and fully mineralized bone, a = 0.7, covered most of the range of values observed in human bone. 
Using dry tissue densities of 1.41 g/cm? for unmineralized osteoid and 2.31 g/cm? for mineralized bone 
(Refs. 16 to 18), a linear approximation for true density through the points (0, 1.41 g/cm’) and (0.7, 
2.31 g/cm’) was obtained: 


p, =1.414+1.290 (5) 


This allowed BVF to be approximated as 


BVF ~ Pee (6) 
Pr 
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Another way to obtain BVF is to piece together the two volume fractions. / can be obtained from the 
following definition. For a given segment or volume of bone, the ash density equals M times the true dry 
density times the maximum ash fraction. 


Pash = M-2.31- (se) (7a) 
= Pash (7b) 
2.31-(Omax) 


Then from the ratio & =P,sp / Papp » O can be solved for from the equation 


—MA2 31 AG) 


a= (8) 
M -2.31+0-1.41 


It is important to note that the value of Omax is assumed to be near 0.7, according to the work of 
Hernandez et al. (Ref. 15), as most of that work, such as Equation (5), is based on the work of Keller 


(Ref. 14). The work of Schileo et al. (Ref. 13), however, revealed the pz p / Papp tatio to be fairly 


constant at a value of 0.6 in human femoral cortical bone. Additionally, their initial analysis showed that 
the average ratio for trabecular bone was lower, 0.46, with a greater variation, 0.34 to 0.62. Consequently, 
Schileo et al. extended their experimental study to include a set of trabecular specimens from the same 
femur that were smaller than those from the initial analysis. For the new specimen set, the average 


Puch b, Papp Tatio was found to be 0.6 with a range of 0.54 to 0.63 and did not differ statistically from that 


of the cortical specimens. Schileo et al. attributed the difference in ratio between the set of smaller 
specimens and the set of initial larger specimens with higher density trabecular bone to possible 
experimental error in the measurement of Papp. There could be an overestimation of Papp in those larger, 
high-density specimen cases, as it is more difficult in the experimental process to assure complete 
removal of water and marrow from the cavities. However, it appeared that the ratio could be accurately 
measured in low-density trabecular specimens, independent of size. In fact, a 0.6 ratio was found for some 
low-density large trabecular specimens. Based on these findings, Schileo et al. rejected their initial 
hypothesis that the ratio decreases as tissue density increases, and a was assumed equal to 0.6 for the 
whole density range.! 


Conversion Methods 


Outlines of the methods are given in terms of initialization of a computational model’s bone volume 
fraction, BVFo value from an initial vVBMDp value. Simulated changes in BVF can then be mapped back to 
vBMD. The first three steps (a, b, and c) are the same for each of the conversion methods. 


Given a VBMD value in terms of pact as measured by QCT, 


a. Convert Pact tO Pash by Equation (2). 
b. Compute Papp by Equation (4). 


c. Calculate initial ash fraction a = Pash F 


p app 


‘A suggestion was made that this could explain the results reported by Hernandez et al. from the large dataset of 
Keller where a low ratio was reported for high-density trabecular bone. 

>There is little difference between Equations (2) and (4). Equation (4) was derived using data from femur 
specimens, which is the site of interest, whereas Equation (2) was derived using data from tibia specimens. 
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Method 1: Uses only the single volume fraction. 


d. Estimate true tissue density p, by Equation (5). 
e. BVF is then given by 


BVF, = Papp (9) 
Pr 


In computational simulations, time-varying changes in BVF; can be converted back to time-varying 
changes in vBMD via changes in ash fraction, apparent density, and ash density. From Equation (9), 
substitute Pash / Papp for ash fraction o in the definition of p; and rearrange to obtain 


BVF\(t)-(1.41-Papp(t) + 1.29- Pash(2)) a (Papp (t))? (10) 


Next, substitute for pasn in terms of Papp from Equation (4) and solve the resulting quadratic for Papp(¢). 
Time-varying changes in ash density Pasn(t) can then be obtained via Equation (4), followed by time- 
varying changes in vBMD(ft) via Equation (2). 

The following methods split BVF into M plus O. 


Method 2: Uses estimate of true density formula, Equations (5) and (6). 


d. Obtain initial volume fraction My from Equation (7b) using OQ 4x =0.7. 
e. Estimate initial true tissue density p, by Equation (5) using initial Oo . 


f. Solve for initial osteoid volume fraction Op) from Equation (6): 


Papp _ Mg-(2.31)+Op-(1.41) _ 


lla 
BVFy Mo +Op me a) 
Oo = 231=Pr -My (11b) 
p, -L.41 
g. BVF is given by 
BVF, =Moy+Op. (12) 


Time-varying changes BV7}(t) = M(t)+O(2) can again be converted back to time-varying changes 
in vVBMD from computational simulations of M(t) and O(f). From Equation (7a), obtain time-varying 
changes in Pasn(¢) and then vBMD(?) via Equation (2). 


Method 3: Uses definition of ash fraction, Equation (8) with On, =0.7. 
d. Obtain initial volume fraction, Mg, from Equation (7b). 
e. Solve for Oo from Equation (8) using o from step c. 


_ Mo: (2.31) -(Gmax- 0) 
(1.41)-ao 


Oo (13) 


f. BVF is given by 
BVF;= My +Op (14) 
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In this case, Mo and Op are provided by the steps under Method 3. Time-varying changes in M(t) and O(2) 
can be used to track a(f) via Equation (8). Then, time-varying changes in vVBMD can be recovered by first 
obtaining Pasn(t) from Equation (7a) and then vBMD(f) via Equation (2). 


Method 4: Uses definition of ash fraction, Equation (8) with O,9, =1.0/1.64. This value of Omax is 
obtained from the maximum value of the ratio 


o. = ———Pash_ (15) 
1.64p asp + 0.01 


as a function of Pash. The steps in Method 4 are the same as those in Method 3, resulting in 
BVF, =My +O (16) 
In this case, Oo is calculated from Equation (13) with the smaller value of Otmax. 


Validation of Bone Volume Fractions Obtained by Conversion 


To test if BVF values calculated by the conversion methods are within appropriate ranges that can be 
expected for healthy adults, the literature was searched for articles that reported BVF values for the 
proximal femur. A study by Tsangari et al. (Ref. 19) reported trabecular BVF values for the 
intertrochanteric region of the proximal femur for male and female adults aged 18 to 49. These results are 
summarized in Table I. 

Two sets of trochanteric vBMD data were subjected to conversion via each method, as summarized in 
Table II. One set represents preflight vBMD data of the trochanteric region in 18 astronauts making ISS 
flights of 4 to 6 months (Ref. 20). The other set represents the pre-bed-rest vBMD of the trochanter region 
in 24 exercisers of a 70-day bed rest study, Countermeasures and Functional Testing in Bed Rest (CFT 
70). QCT vBMD data was provided by the NASA Johnson Space Center Bone Lab through the NASA 
Life Sciences Data Archive. 

Comparing the calculated BVF values to the experimental values presented in Table I, Methods 2 and 
4 present an average and standard deviation for the CFT 70 bed rest group that is inside and toward the 
lower half of the range for both the male and female experimental averages and standard deviations. For 
the astronaut group, calculated BVF values for Methods 2 and 4 are close to the female experimental 
values in terms of average and standard deviation and are fairly centered inside the experimental range of 
values for males. The calculated values of Methods | and 3 for the astronaut group lay inside of the male 
experimental range in terms of average and standard deviation and are shifted toward the higher half. 


TABLE I—EXPERIMENTAL VALUES FOR TRABECULAR 
BONE VOLUME FRACTION (BVF) FOR PROXIMAL 


FEMUR REPORTED IN REFERENCE 19 
Gender Age range Trabecular BVF 
Female n= 13 18 to 49 0.146+0.04 
Male n= 14 18 to 49 0.151+0.064 
TABLE II.—PREDICTED TRABECULAR BONE VOLUME FRACTION (BVF) VALUES FOR EACH METHOD 
Trabecular VBMD Age range Calculated BVF values 
data source Method 1 Method 2 Method 3 Method 4 
a 
aA eae Bue 0.145+0.02 0.136+0.02 0.149+0.02 0.139+0.02 
n=24 3245.3 
ae ie 9) 44.6+4.0 0.151+0.017 0.142+0.016 0.156+0.017 0.145+0.016 


*Countermeasures and functional testing 
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Since the standard deviations of the calculated BVF values for each method are small, a clearer 
comparison between the calculated BVF values of the methods and the experimental BVF values can be 
made by observing the distribution of the individual BVF values of the methods. These are illustrated in 
Figure | for the CFT 70 bed rest subjects and in Figure 2 for the astronauts. 

This makes a much clearer evaluation of the methods. As shown in Figures | and 2, the distribution 
of the individual BVF values falls within the male experimental range for all four methods. For the bed 
rest subjects, distributions of Methods 2 and 4 are shifted toward the lower end of the experimental 
ranges, with several values outside the lower end of the female experimental range, which explains the 
lower averages in Table II. However, as the bed rest subjects are all males, their BVF values are all within 
the experimental range for the four methods. 


+ BVF, 
4 * * RR EK RR KK KKK * peal: 
© BVF3 
* BVF4 
3} © © oO 20 OO KY 00 9 900 © ; 


2 0 Oo O CM@® O CO CD COO 000 O 4 
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ee ties 
Female BVF range 
Male BVF range 


Figure 1.—Distribution of bone volume fraction (BVF) values of individual bed rest subjects for 
Methods 1 to 4. The range of experimental values for females and males according to the 
average and standard deviation are shown below the horizontal scale of BVF values. 
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Figure 2.—Distribution of bone volume fraction (BVF) values of individual astronauts for 
Methods 1 to 4. The range of experimental values for females and males according to the 
average and standard deviation are shown below the horizontal scale of BVF values. 
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TABLE IIT —MINERALIZED AND OSTEOID VOLUME FRACTIONS FOR 
AVERAGE ASTRONAUT vBMD = 0.148 g/cm? 


Method 2 Method 3 Method 4 omy =1/1.66 Method 4 an, = 0.60103 
M O M O M O M O 
0.1203 | 0.0217 | 0.1203 | 0.0363 | 0.1381 0.0071 0.1401 0.0038 


Compared with the distribution of BVF values of the bed rest subjects, the distribution of the BVF 
values of the astronauts is shifted toward the higher end of the experimental ranges for each method. This 
is reasonable, as the astronauts are a healthier group that undergoes strength conditioning and training 
prior to flight. 

The same article that reports intertrochanteric BVF values (Ref. 19) also has mean O values. 
Compared with M, this normally tends to be small. The reported values are on the order of 10% for 
females between 18 and 49 and 10° for males between 18 and 49. Table III shows the values of Methods 
2 to 4 for the average astronaut vBMD of 0.148 g/cm’. 

For the last set, the value of 0.60103 is obtained from Equation (15) by setting the value of Pash to 
1.2, which was usually the estimated top value of Pash in Reference 12. 

For consistency, modeling that uses the sum of mineralized and osteoid volume fractions should be 
the same or close to the values of the single volume fraction. To achieve that, Methods 2 and 4 can be 
adjusted. 

For Method 2, consider the following alternative: 


Method 2a: 


d. Obtain initial volume fraction, Mo, from Equation (7b). 
e. Set Mo+Oo = BVF;, defined in Equation (9). 
f. Solve for Oo. 


Op = PHP — Mo (17) 
Pr 


BVFya =Mo +Opo 


g. (18) 


BVF»5a will then have the same values as BVF}. 

For Method 4, the results of the study by Schileo et al. (Ref. 13) for smaller trabecular specimens can 
be considered. In those cases, some variation in ash fraction up to 0.63 was observed. This suggests using 
an estimate of Omax of about 0.63 to 0.65. So the following alternative for Method 4 is defined: 


Method 4a: 


d. Obtain initial volume fraction, Mo, from Equation (7b) with Omax = 0.65. 


Pash 
M = 19 
2.31-(0.65) am” 


e. Solve for Oo from Equation (8). 


— Mo -(2.31)-(0.65-a9) 
7 (1.41)- a 


(20) 


f. BVFya = Moy + Op (21) 
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For Method 4a, BVFya =0.144+0.02 for the bed rest subjects, which is similar to BVF in Figure 1. 
For the astronauts, BVFya =0.150+0.019, which is similar to BVF; in Figure 2. 

Alternately, Method 1 can be adjusted. Since the work of Keyak et al. and Schelio et al. suggests a 
relatively constant ash fraction value of 0.6, let 1/1.64 = 0.6097 be an estimate of the maximum value of 
Omax from Equation (15) and replace Equation (5) with a line through (0, 1.41g/cm*) and (0.6097, 
2.31g/cm’). This will result in an equation that will give a larger value of true density p, in Equation (9) 
and hence a smaller value of BVF;. This will generate values of BVF; close to the values of Method 4. 


Bone Volume Fraction From Microcomputed Tomography 


A study by Tassani et al. (Ref. 21) to investigate whether tissue mineral density can be assumed a 
constant in adult human bone used microcomputed tomography (mCT) analysis. Ninety-six specimens 
from various lower limb sites (tibias and femurs) were extracted from two female cadavers. BVF was 
measured via mCT, and specimens were ashed to gravimetrically measure Pash. Performance of a 
regression analysis between Pash and BVF for the pooled groups of trabecular and cortical specimens 
obtained the following relation: 


Pash =1.18BVFinc, +0.01 (22a) 


An intrasite analysis was performed on a second group of 19 trabecular specimens extracted from 
femoral heads of different donors. Performance of a regression analysis of Pasn over BVF gave the relation 


Pash = 0.99 BVFinc, + 0.04 (22b) 


If the data were to be used to obtain BVF values directly from Pash, one might solve for BVFinct in 
either equation to get a valid estimate. However, forming a regression of BVF on Pash would not generally 
give the same equation. Unfortunately, the data was not available to perform the regression analysis that 
way. If Equation (22a) is reversed, values of Pash from Equation (2) generate values of BVFinct = 
0.156+0.019 for the astronaut’s vBMD and BVFnct = 0.149+0.02 for the bed rest subject’s vBMD. Values 
generated from the reverse of Equation (22b) and the reverse of the regression from the trabecular group 
of the original 96 specimens produced very similar results. These are average values similar to those 
values listed under Method 3 in Table II. It would be worthwhile to obtain the data from the study by 
Tassani et al. and perform a regression of BVF on Pash to obtain the correct reverse of Equations (22a) or 
(22b), then substitute the values of Pasn from Equation (2) into the reverse of Equation (22a) or (22b) and 
compare the results to the other methods. This would be the simplest and most direct method of 
conversion. 


DAP Computational Bone Physiology Method 


Currently the DAP model that predicts changes in BVF uses conversion Method 2 for converting 
vBMD into BVF by combining M and O. This early version of the conversion is based on knowledge 
available at the time of the work of Keyak et al. (Ref. 12) and Hernandez et al. (Ref. 15). This method 
was used in the validation analysis of the femoral neck model. At that time, the intention was to use a 
modulus formula developed by Hernandez et al. in the finite element model for bone that is coupled with 
the model’s bone physiology component. In the process of validating the model for the femoral neck, the 
modulus formula was changed to a formula by Keyak et al. (Ref. 11). Further literature research yielded 
additional knowledge of the other cited work. This, along with continuing mathematical research into 
other conversion possibilities, and prompted by a desire to extend the model to the full proximal femur, 
led to consideration of the other methods. As a result, researchers plan to modify the code to use Method 
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4 for consistency, as it makes use of the Keyak and Schileo ash fraction values of about 0.6. Method 4 
appears to validate O values better, as indicated by Table HI. Also planned is completion of an updated, 
shorter simulation in which BVF changes are not split into M and O. 


Conclusions 


The methods outlined here for converting bone mineral density (BMD) derived via quantitative 
computed tomography (QCT) to an approximated bone volume fraction (BVF) in trabecular bone of the 
proximal femur are based on multiple experimental studies, correlation relations derived from those 
studies, and definitions of bone properties. The best way to validate the methods would be to apply the 
conversion methods to the BMD values of the subjects used to obtain the experimental ranges of BVF 
values, but that data is unavailable. Although the BVF values predicted by the methods differ, they all fall 
within a range of experimental values for the intertrochanteric region. Judgment of the methods was thus 
limited to how well the distribution of BVF values obtained from pre-bed-rest BMD of 24 subjects and 
preflight BMD of 18 astronauts covers the experimental range. The need for Digital Astronaut Project 
model changes in BVF derived by Method 1 to be consistent with changes in BVF derived by methods 
that combine M and O, such as Methods 2, 3, and 4, was also considered. Ultimately, with any method, 
the interest will be in validation of the computational model. That involves an analysis of how well the 
computational model’s resulting changes in VBMD compare with experimental changes in VBMD. 
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Appendix—Definitions 


e Apparent density—Weight of bone tissue divided by total volume of segment or specimen 


e Real density—Weight of bone tissue divided by total volume of bone tissue (real volume) 
© Can be referred to as true density 
o Compact bone approximately constant—1.9 kg/cm? (Ref. 22) 


e Volumetric bone mineral density (VBMD)—Bone mineral content divided by total volume of bone 
segment. Measure of bone density reflecting strength of bones as represented by calcium content 


e Ash density (apparent)—Ratio between ash weight and volume of sample of bone segment or 
specimen 


e Bone volume fraction (BVF)—Volume of bone tissue divided by total volume of bone segment 


© Bone volume includes internal pores like lacunae and canaliculi. 
o Trabecular bone may have BVF ranging from just over 5 percent to a maximum of 60 percent. 


e Ash fraction—Ashed mass of a bone sample divided by the dry mass of the same sample 

© Sometimes referred to as percent mineralization 
Considered a standard for determining amount of mineralization 
Has been used to estimate bone material properties such as elastic modulus (Ref. 14) 
Calculation requires destructively ashing the specimen 
The mathematical expression of the ratio can vary, although the expressions have the same 
meaning: 


OoO00 0 


ash fraction = ash mass/dry mass 


= inorganic mass (mass of Ca hydroxyapatite) 


inorganic mass + organic mass 


= mineralized mass 


mineralized mass + unmineralized mass 
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